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Abstract 

We discuss the solution of a nonlinear ordinary differential equation that appears in 
a model for MHD viscous flow caused by a shrinking sheet. We propose an accurate 
numerical solution and derive simple analytical expressions. Our results suggest 
that a recent perturbation treatment of the same problem exhibits a pathological 
behaviour and conjecture its probable cause. 



1 Introduction 

There has recently been enormous interest in the so called homotopy per- 
turbation methods. In particular there is an open controversy between the 
developers of homotopy perturbation method (HPM) and homotopy analysis 
method (HAM)[1,2J. One of the controversial points is related to the appear- 
ance of some "secular" terms in the perturbation solutions. Whereas He[T] 
proposed their removal, Liao[2] argued that they are harmless because they 
vanish as the variable increases towards infinity. 
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Sajit and Hayat[3] have recently discussed the MHD viscous flow due to a 
shrinking sheet. They converted the model partial differential equations into a 
nonlinear ordinary differential equation which they solved by means of HAM. 
The authors discussed the convergence of the perturbation approach and ob- 
tained apparently accurate results for several values of the model parameters. 
In fact, they concluded that "The obtained HAM solution is valid for all values 
of the suction parameter and Hartman number." 

The purpose of this paper is to analyze the accuracy of those HAM results 
because we believe that they may cast light on the abovementioned discussion 
about the "secular" terms. In Sec. [2] we consider an accurate numerical calcu- 
lation, in Sec. [3] we derive simple approximate analytical expressions, and in 
Sec. S] we discuss the results for a particular choice of the model parameters 
and draw conclusions. 

2 Accurate numerical calculation 

Since we are not interested in the validity and usefulness of the model we 
just concentrate on the equation that the authors solved approximately by 
means of HAM. By means of an appropriate transformation Sajit and Hayat[3j 
converted the model partial differential equations into the ordinary nonlinear 
differential equation 

f"'( V ) - M 2 f'( V ) - f'( V ) 2 + mf(r))f"(r)) = 
/(0) = s, /'(0) = -1 

\vmf\ V )=0 (1) 

where a prime indicates differentiation with respect to the variable rj and M, 
m and s are model parameters. The main problem is to obtain the value of 
f"(0) that is consistent with the condition at infinity. Once we have it, then 
we can resort to any numerical integration routine to obtain the solution for 
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all rj > 0. 



In order to determine f"(0) = a accurately we apply the Hankel-Pade method 
developed some time ago that proved successful for the treatment of two- 
point boundary value problems [i|5|6j . It consists of expanding the solution in 
a Taylor series about 77 = 



and then calculating the roots of the Hankel determinant H D (p) with ma- 
trix elements fi + j + d(ct), i,j = 1,2, ... ,D. The calculation is straightforward 
because the matrix elements and, consequently, the Hankel determinant are 
polynomial functions of the unknown parameter a. As shown in earlier ap- 
plications of the method[lj5 ; 6j one expects to find a sequence of roots an, 
D = 2,3, . . ., that converges towards the appropriate value of /"(0). As said 
above, once we have a sufficiently accurate value of this parameter then we 
can apply any numerical integration algorithm and obtain f(rj) for all rj > 0. 
Alternatively, in some cases the Pade approximants, on which the method is 
based, give sufficiently accurate results[5]. 

3 Approximate analytical expressions 

Approximate analytical solutions to the equations of a physical model com- 
monly provide greater insight into the nature of the phenomenon under in- 
vestigation. The analytical expressions provided by HAM[3J appear to be so 
complicated that one can only use them within a computer algebra system. 
In this sense this kind of solution is not much different from the results pro- 
vided by the numerical integration routines that are also built in most such 
software packages. The purpose of this section is to provide simple analytical 
expressions for the straightforward discussion of the MHD model. 



00 




(2) 



j=0 



3 



As 7] increases fit]) 2 is expected to be smaller than the other terms and 
therefore f(rj) will behave approximately as the solution of f'"{rf) — M 2 f'(r)) + 
m foof"(v) = 0; where = f(i] — > oo). In other words, we expect that 
f'(rj) ~ be~~P v for sufficiently large rj. Therefore, it seems reasonable to try the 
ansatz 

N 

fW(ri) = ^bje-**, N = 1,2,... (3) 

j=0 

If we substitute it into the differential equation (pQ) we obtain an expression 
of the form 

27V 

Y / Rj(b 1 ,b 2 ,...,b N ,P)e-^ = (4) 
i=i 

Therefore, the optimal values of the adjustable parameters bj and (3 should be 
solutions to 

/M(0) = «, fW(o) = -i 

Rj(b u b2,---,bN,P)=0, j = l,2,...,iV (5) 



Since it is our purpose, for the reasons already given above, to keep the re- 
sults as simple as possible we just consider the first two approximation orders 
explicitly. For N = 1 and N = 2 we easily obtain 

1 , 1 n y/4M 2 + m?s 2 -Am + ms) . . 

b = s--, bl = -^= (6) 



and 



f3 2 -M 2 2(M 2 - f3 2 ) + m(2(3s - 1) /3 2 — M 2 + m(l — /3s) 

00 - n > °1 - o ' ° 2 _ o 

mp mp mp 

4/3 4 - 4/3 3 ms(2 - m) - 2 / 3 2 [2m(mV - ms 2 - 1) - M 2 (3m - 4)] 
-2f3ms[M 2 (5m - 4) - 2m(m - 1)] 

-2M 4 (3m - 2) - 2M 2 m(2 - 3m) - m 2 (m - 1) = (7) 



respectively. 
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It is probable that this approach, or somewhat similar to it, had been used in 
the past. However, it is useful for our purposes and we are not aware that it 
had been applied to the problem discussed by Sajid and Hayat[3]. 

4 Results and discussion 

Sajid and Hayat[3] analyzed the form of f'{rf) for several values of the model 
parameters. Here we simply consider the case M = m = 2 for the largest 
values of s for which f'(rj) appears to exhibit a maximum according their 
Figure 2. 

For s — 1.8 a sequence of roots of the Hankel determinant Hp(a) with 
d = 1 and D = 2, ... ,30 suggests that a = 4.20411340 (the reader may find 
examples of the rate of convergence of the method elsewhere [4"|5|6] ) . Sajid 
and Hayat[3j only showed values of f"(0) for M = 2, s = 1 and m = 1,2; 
therefore, their result for s = 1.8 is not available for comparison. The correct 
behaviour of the numerical Runge-Kutta solution for large values of r\ pro- 
vides an additional confirmation of the accuracy of that value of /"(0). Fig. [1] 
shows such numerical results and also those given by the analytical expressions 
([6]) and (TH). Our simple analytical expressions already provide satisfactory re- 
sults for all values of r\ as well as the following acceptable estimates of /"(0): 
a' 1 ' = V131/5 + 9/5 w 4.1 and cJ 2 ' w 4.198. We appreciate that the accuracy 
increases with N (at least for the first two approximations). Besides, the fact 
that bi ~ 0.238 ^> b 2 ~ 0.00309 for N = 2 suggests a remarkable convergence 
rate. 

Fig. Q] clearly shows that neither the approximate analytical expressions nor 
the accurate numerical results exhibit a maximum. Therefore, we conclude 
that the maxima found by Sajid and Hayat[3] when increasing the suction 
parameter are merely artifacts of the HAM. In our opinion, such spurious 
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maxima are probably caused by the "secular" terms of the form r] k e~ nv shown 
in Eq. (24) of Ref. [3] and discussed by He an Liao fTp] . Although those terms 
certainly vanish as 77 — > 00 they may have some non-negligible undesirable 
effect for moderate values of r\. 

Summarizing: we have verified once more that the Hankel-Pade method is 
useful for the treatment of two-point boundary value problems by obtaining an 
accurate value of the unknown parameter appearing in the nonlinear equation 
for a MHD viscous flow model[3]. It is necessary for a successful application 
of any numerical integration routine. We have also derived simple accurate 
analytical expressions that may be useful for the discussion of the physics of 
the problem. Both the analytical expressions and the numerical results have 
revealed a pathological behaviour of the results produced by the much more 
elaborate approach called HAM[3J. In this way we hope to have settled the 
argument about the "secular" terms in the HAM expressions [TP] . Perhaps, 
the HAM users may want to verify our conjecture and throw some more light 
on the subject. 
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Fig. 1. Numerical results (circles), first-order (dashed line) and second-order (solid 
line) analytical expressions for f'(r)). 
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